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ABSTRACT 


The  objective  of  the  research  has  been  to  develop 
engineering  methodologies  applicable,  but  not  limited,  to 
aerospace  automatic  control  design  problems  in  which  there  are 
performance  specifications  requiring  precise  control  of  system 
behavior  in  the  presence  of  stochastic  disturbances  (e.g.,  wind 
gusts)  and  large-but-bounded  uncertainties  in  the  dynamical 
response  of  the  system  (e.g.,  parameter  uncertainty,  unmodeled 
nonlinearities,  and  so  forth)  During  the  past  three  years  of 


AFOSR  supported  research,  fk  cohesive  body  of  theory  has  been 
developed  that  enables  engineers  to  relate  the  ability  of  feed¬ 
back  control  systems  to  meet  such  specifications  directly  and 

S'  9^ 

quantitatively  to  the  ^return  difference  matrix**  associated  with 

ft  in-Cftity 

the  system's  feedback  loops.  Now  results  enablings'll  optimi- 
ration*^  of  returned  difference  singular  value  Bode  plots  pro¬ 
mise  to  be  of  great  value  in  robust  multivariable  feedback 


controller  synthesis.  Continuing  research  is  currently  being 
aimed  at  further  tightening  the  links  between  this  theory  and 
the  most  recent  developments  of  modern  stochastic  linear  optimal 
control  synthesis  thoery,  and  extending  the  results  to  admit 
more  practical  problems,  so  that  this  thoery  may  be  used  more 
effectively  by  engineers  to  efficiently  and  systematically  design 
the  feedback  gains  that  determine  a  feedback  system's  return 
difference  matrix.  Such  results  substantially  reduce  the  depen¬ 
dence  of  control  engineers  on  intuition,  simulation,  and  luck 
and  provide  the  know-how  to  successfully  and  efficiently  solve 


INTRODUCTION 


The  increased  demands  for  quick  and  precise  control  over 
aircraft  and  space  vehicle  response  that  are  anticipated  in 
the  coming  decades  will  have  to  be  matched  with  automatic  con¬ 
trol  systems  that  can  respond  instantaneously,  without  moment- 
to-moment  human  guidance,  anticipating  vehicle  response  inso¬ 
far  as  is  possible  and,  more  importantly,  continually  and  auto¬ 
matically  monitor  vehicle  response  and  re-adjust  control  sig¬ 
nals  to  correct  for  unpredicted  deviations  from  the  desired 
response.  Such  unpredicted  response  variations  can  result  from 
external  disturbances  (e.g.,  wind  gusts)  and  from  the  impossibi¬ 
lity  of  employing  a  sufficiently  complex  and  accurate  model  of 
the  vehicle's  dynamics  to  account  for  every  vibrational  mode, 
every  nonlinearity,  ...,  every  variable  affecting  system  res¬ 
ponse  . 

While  the  state-space-based  mathematical  theory  for  control¬ 
ling  systems  without  substantial  uncertainty  regarding  dynamical 
response  grew  relatively  sophisticated  during  the  two  decades 
of  the  1960's  and  1970's,  there  was  almost  no  significant  prog¬ 
ress  concerning  the  control  of  systems  having  uncertain  response 
since  the  1940 's  and  1950's  when  great  strides  were  made  in  the 
development  of  the  "classical"  transfer-function-based  theory 
for  the  control  of  simple  single-input-single-output  (SISO) 
linear  time-invariant  (LTI )  systems  with  uncertainty.  Conse¬ 
quently,  when  this  research  project  was  begun  in  October  1979 
there  was  no  adequate  theory  to  provide  engineers  with  an  effi- 
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cient,  systematic  procedure  for  the  design  of  precision  con¬ 
trollers  for  more  complex  multi-input-multi-output  systems 
such  as  the  highly  unstable,  fast  responding,  control  confi- 
fured  aerospace  vehicles  that  are  expected  to  be  operating 
in  tomorrow's  combat  environment.  The  objective  of  the  pre¬ 
sent  research  has  been  to  develop  this  badly  needed  theory 
so  that  the  engineers  who  must  design  tomorrow's  aerospace 
vehicles  will  have  more  than  intuition,  trial-and-error  simula¬ 
tion  and  mystical  "seat-of-the-pants"  insight  for  guidance 
in  designing  uncertainty-tolerant  automatic  controllers  for 
these  vehicles. 


5 


BACKGROUND  AND  PROGRESS 


Since  work  began  on  this  project  in  October  1979,  the  re¬ 
search  effort  has  been  generally  successful  in  achieving  its 
main  objective  of  relating  the  return  difference  matrix  to 
the  uncertainty  tolerance  properties  of  a  system.  These  pro¬ 
perties  are  also  known  as  robustness  properties  or  feedback 
properties.  Our  results,  together  with  some  related  results 
useful  in  the  actual  synthesis  of  robustly  uncertainty  tole¬ 
rant  feedback  controllers,  were  reported  in  the  paper  "Feed¬ 
back  Properties  of  Multivariable  Systems:  The  Role  and  Use 
of  the  Return  Difference  Matrix"  [1  ]  .  This  paper  discusses 
the  central  roles  in  feedback  theory  of  the  return  difference 
matrix  (denoted  I+L(s))  and  the  inverse-return  difference  mat¬ 
rix  (denoted  I+L-1 (s) ) . 

Among  the  new  theoretical  results  in  [ 1  ]  are  the  following: 

(i)  A  new  method  for  exact  evaluation  of  the  sensitivity 
of  multivariable  feedback  control  systems  which  overcomes  sig¬ 
nificant  practical  limitations  associated  with  previously  known 
methods.  Sensitivity  to  large  plant  and  sensor  variations  is  directly 
related  to  the  nominal  system's  return  and  inverse-return  difference  matri¬ 
ces,  (See  Theorem  2.1  and  Theorem  2.2  in  [  1  ]  )  . 

(ii)  Significant  drawbacks  of  characteristic  locus  analy¬ 
sis  methods  (cf.  MacFarlane  and  others)  are  described.  Return 
and  inverse-return  difference  singular  value  plots  are  found  to 
overcome  some  of  the  drawbacks  of  characteristics  loci.  The 
results  have  been  found  to  be  useful  in  quantifying  some  funda¬ 
mental  limits  on  the  achievable  performance  of  feedback  control 
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systems.  (See  Section  3  and  Section  4  of  [  1  ]  ) . 

(iii)  A  technique,  based  on  stochastic  linear  quadratic 
Gaussian  (LQG)  optimal  control  theory,  has  been  developed  to 
aid  the  shaping  of  the  return  and  inverse-return  difference 
singular  value  plots.  Though  the  technique  is  to  a  certain 
extent  a  trial-and-error  design  technique,  it  continues  to  be 
substantially  more  systematic  than  any  other  method  that  is 
currently  available  for  synthesizing  multivariable  control 
systems  to  meet  specifications  requiring  a  robust  tolerance 
of  disturbances,  noise  and  plant/sensor  modeling  errors. 

(See  Section  5  in  [  1  ]. )  To  demonstrate  the  viability  of  the 

of  the  technique,  the  theory  has  been  applied  to  the  synthe¬ 
sis  of  an  automatic  controller  for  the  longitudinal  dynamics 
of  an  advanced  control  configured  vehicle  (CCV)  aircraft,  viz., 
the  NASA  HIMAT  remotely  piloted  aircraft  (see  Section  6  of 

fl  ]  )  . 

More  recently  research  effort  has  focused  on  several  issues 

First,  substantial  effort  has  been  focused  on  the  important 
practical  issue  of  how  to  solve  the  so-called  "Inverse  Problem 
of  Linear  Quadratic  Gaussian  (LQG)  Optimal  Control"  in  the  aene 
ral  setting  in  which  the  controller  is  dynamical  and  the  plant 
is  subject  to  plant  and  sensor  noise.  This  inverse  problem  is 
as  follows:  Given  a  realizable  closed-loop  control  return-dif¬ 
ference  matrix,  the  plant  transfer  function  matrix  P(s),  and 
.!  e  •  ant  and  sensor  noise  power  spectra  matrices,  say  <  (s) 

and  In(s) ,  find  the  linear  quadratic  cost  matrices  F(s)  and 
Q(s)  such  that  the  closed-loop  system  is  optimal  in  the  sense 
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that  the  following  stochastic  cost  is  minimized; 


■  -f 

2”J  J 


Tr(t!(s)r  (s)  +  R(s)  Iu(s))  ds 


where  Ey(s)  and  £u(s)  denote  the  respective  closed-loop 
power  spectra  matrices  of  the  sensor  output  and  control  in¬ 
put  signals.  We  are  pleased  to  report  that  the  solution  to 
the  problem,  together  with  extensive  discussion  of  its  ra¬ 
mifications  regarding  a  number  of  related  problems  is  in 
hand  [5 3  .  This  work  was  supported  under  AFOSR  Grant  80-0013. 

The  importance  of  the  LQG  inverse  problem  to  the  proper 
understanding  of  robust  multivariable  feedback  control  system 
design  cannot  be  understated.  LQG  multivariable  feedback  de¬ 
signs  are  preferred  for  a  variety  of  reasons:  Good  computer 
software  is  readily  available  for  solving  the  LQG  design  equa¬ 
tions  (e.g.,  [10]  ) ,  LQG  designs  optimize  inherent  trade-offs 
between  robustness  properties  such  as  sensitivity  versus  sta¬ 
bility  margin  [l]  ,  and  many  engineers  in  the  aerospace  indus¬ 
try  are  familiar  with  the  basic  LQG  concepts.  The  solution  to 
the  inverse  problem  plays  a  vital  role  in  our  understanding 
of  how  one  can  use  the  LQG  theory  to  place  the  poles  and  zeros 
of  3>L(s)  in  order  to  achieve  acceptable  robustness  and  sensi¬ 
tivity  singular-value  Bode  plots  [l]  and  acceptable  transient 
response  and  asymptotic  tracking  properties  [5  ]  .  Also,  as  it 
is  common  practice  to  iteratively  adjust  the  LQG  cost  matrices 
when  "fine-tuning"  a  control  desian  to  meet  various  robustness 


specifications,  the  solution  to  the  inverse  problem  provides 
a  starting  point  for  fine-tuning  non-LQG  multivariable 
feedback  designs. 

A  second  focus  of  the  research  effort  is  one  of  rather  funda 
mental  significance  in  linear  feedback  theory:  realizability. 

A  closed-loop  feedback  system's  return-difference  I  +  L(s)  is 
said  to  be  realizable  for  a  given  plant  P(s)  if  for  some  con¬ 
troller  C(s) 

(i)  L (s)  =  P(s)C(s)  ,  and 

(ii)  (I  +  L(s))-1  is  stable,  and 

(iii)  P  (s)  and  C(s)  are  "coprime"  in  the  sense  that  there 
are  no  unstable  pole-zero  cancellations  in  the  product  P(s)C(s). 

A  new  simplified  realizability  result  (Lemma  1,  [4  ]  )  has  been 

developed  which  characterizes  realizability  directly  in  terms 
of  the  poles  and  associated  residues  of  the  return-difference 
and  inverse-return  difference.  This  is  an  improvement  over 
previous  multivariable  results([  3  ],  Lemma  3)  which  require 

a  solution  of  the  so-called  "Bezout"  equation  and  cive  only 
limited  insight  into  the  constraints  imposed  by  realizability 
on  the  set  of  achievable  return  difference  anc  inverse-return 
difference  matrices. 

The  realizability  question  arose  in  connection  with  the 
previously  mentioned  LOG  inverse  problem,  and  our  new  reali¬ 
zability  result  in  [4]  plays  a  crucial  role  in  the  solution 
of  the  LOG  inverse  problem  in  [5  ],  in  addition  to  layinc 
the  groundwork  for  stability  marain  optimization  problem  solved 
in  [4  ] . 


The  problem  of  stability  margin  optimization  and  the 
mathematically  dual  problem  of  sensitivity  minimzation  with 
respect  to  the  L  norm  have  been  a  third  major  thrust  of  our 
research  effort  in  the  past  year.  With  the  aid  of  results 
in  [l  ]  which  show  that  stability  margin  is  inversely  pro¬ 
portional  to  the  L°°  norm  of  the  inverse  of  the  inverse-return 
difference  matrix  I  +  L~x (s)  and  the  aid  of  our  new  reali¬ 
zability  result  (Lemma  1  of  [4  ]  ) ,  it  is  shown  in  [4]  that 

the  problem  of  designing  a  feedback  controller  to  maximize 
stability  margin  (subject  to  decoupling  and,  perhaps,  asymp¬ 
totic  tracking  constraints)  is  mathematically  equivalent  to 
the  minimization 

min  lix  •  (s)  lloo  ,  (j  ■  1,  - -  m) 

Xj  (s)  €  H°° 

subject  to  complex  interpolation  constraints  of  the  for™ 

Xj(sij)  =  wij,  (i  =  1,  ...,  n)  ,  Refs^)  >  0, 

where  s^j  and  w^j  are  complex  constants  and  K00  is 

the  Hardy  space  of  stable  transfer  functions  with  the  I & 
norm , 

!.*(s)  L  £  sup|x(  jw)  I . 

A  simple  solution  to  this  interpolation  problem  requir.nq 
only  the  calculation  of  certain  eigenvectors  and  eigenvalues 
is  available  in  the  mathematics  literature  [39  ] .  This  leads 
in  turn  to  the  multivariable  feedback  controller  havina  maxi¬ 
mal  stability  robustness.  It  also  points  the  way  to  imrravinc 


and  extending  to  multivariable  systems  certain  recent  results 
of  Zames  and  Frances  [40]  concerning  single- loop  feedback  sen¬ 
sitivity  minimization  with  respect  to  the  L“  norm.  The  results 
recently  have  been  generalized  to  the  decoupled  multivariable 
case  by  the  principle  investigator  and  Ph.D  student  B.S.  Chen  [4,5]. 

Another  area  which  we  have  studied  during  the  past  year 
has  been  the  use  of  the  so-called  Hankel  singular  values  [8-15] 
as  an  alternative  measure  of  the  "sizes"  of  the  return  difference 
I  +  L  (s)  and  inverse  return  difference  I  +  L~*(s)  for  purposes  of 
robustness  continuation  proposal,  the  Hankel  singular  values 
are  attractive  for  at  least  two  reasons: 

(1)  There  are  finitely  many  Hankel  singular  values  for  a 
finite  order  system  whereas  conventional  frequency  response 
singular  values  involve  an  infinite  continuum  of  w-dependent 
singular  values  (several  at  each  frequency  a  €  (-  00 ,  °°)) .  It  can 
be  shown  that  this  discrepancy  stems  from  the  following  funda¬ 
mental  mathematical  fact  -  the  Hankel  operator  is  compact  with 
a  finite,  discrete  spectrum  (i.e.,  set  of  eigenvalues),  whereas 
the  frequency-response  method  corresponds  to  noncompact  Toeplitz 
operator  with  a  continuous  spectrum  (with  eigenvalues  equal  to 
the  set  eigenvalues  of  the  operators  frequency  response  matrix 
as  frequency  o>  ranges  over  all  real  values  u  €  [- <*-  ,’*)). 

(2)  Hankel  singular  value  decomposition  of  a  linear  sys¬ 
tem  captures  the  causal  nature  of  a  system  and  is  closely  re¬ 
lated  to  well-understood  state-space  concepts  such  a  controll¬ 
ability  and  observability.  In  contrast  frequency  response 
singular  values  do  not  enable  one  to  readily  distinguish  realiz¬ 
able  (i.e.,  causal)  from  non-realizable  (i.e.,  non-causal  systems 


On  the  encouraging  side,  Hankel  singular  values  have  been 
found  to  correlate  well  with  the  practically  relevant  time- 
response  and  frequency-response  properties  of  lightly  dam¬ 
ped  aerospace  structures  [14  ]  .  On  the  other  hand,  we  have 

found  that  one  can  synthesize  classroom  examples  of  transfer 
function  pertubations  which  have  vanishingly  small  Hankel  sin¬ 
gular  values,  yet  are  capable  of  destabilizing  the  most  ro¬ 
bust  feedback  design;  i.e.,  Hankel  singular  values  seem  to 
induce  an  inappropriate  topology  in  the  space  of  transfer  mat¬ 
rices  [15  ]  . 

Still  there  are  close  links  between  Hankel  singular  value 
optimization  problems  and  L00  norm  frequency  response  robust¬ 
ness  optimization  problems  involving  the  return  difference  and 
inverse-return  difference.  The  stability-constrained  L00  op¬ 
timal  complex  interpolation  problems  that  yield  to  the  solution 
of  sensitivity  minimization  problems  [40]  and  stability  rnarcin 
maximization  problems  [39]  also  provide  the  solution  to  rela¬ 
ted  Hankel  singular  value  optimization  problems  involving  re¬ 
duced  order  system  modeling  (  Pee  [  9]  and  [  41  ]  ).  Thus, 
although  we  do  not  presently  see  any  future  for  the  Hankel  sin¬ 
gular  values  as  a  measure  of  robustness,  it  is  clear  now  that 
concepts  and  theoretical  results  developed  in  conjunction  with 
optimal  Hankel  model  reduction  theory  will  be  useful  to  us  in 
performing  the  realizability-constrained  L00  -  norm  frecuency- 
response  singular  value  optimizations  that  naturally  arise  in 
synthesizing  optimally  robust  multivariable  feedback  controDler 
designs. 


CURRENT  RESEARCH 


The  primary  focus  of  our  current  research  is  on  the  synthesis 

of  optimally  robustly  uncertainty-tolerant  multivariable  feedback 

control  systems.  The  results  presently  in  hand  for  this  purpose 
2 

include  the  L  -norm  LQG  optimal  synthesis  [1,  Section  5],  L°°-norm 
optimal  sensitivity  synthesis  [40]/  and  L°°-norm  optimal 

stability  margin  synthesis  [4].  The  LQG  approach  developed 
in  [l]  is  at  present  the  most  useful  of  the  three  because  it 
allows  one  to  establish  an  optimal  trade-off  between  the  us¬ 
ually  conflicting  robustness  requirements  of  sensitivity  ver¬ 
sus  stability  margin.  One's  choice  of  LQG  cost/noise  weighting 
matrices  allows  one  to  specify  in  each  instance  how  much  im¬ 
portance  is  to  be  attached  to  each  of  these  two  conflicting  de¬ 
sign  requirements.  However  the  Is -norm  which  is  the  basis  for 
LQG  optimization  is  inappropriate  for  stability  margin  analy¬ 
sis  since  even  plant  perturbations  Ap(s)  having  vanishingly 
small  L*-norm, 


Tr  (b 


<3^  )  L(  j  -  )  ) 


can  cause  instability;  i.e,  ,  the  L^,  norm  like  the  Hankc-l  —norm 
[15  ]/  gives  an  inappropriate  topology  to  the  space  or  transfer 


matrices  for  stability  margin  analysis.  For  stability  margin 
analysis  the  L°°  norm  associated  with  frequency  response  sin¬ 
gular  values  is  known  to  be  much  more  appropriate  [  1,  17,  42, 
43  ].  For  sensitivity  optimization  the  advantages  of  the  L°° 
norm  over  the  La  norm  are  not  so  clear-cut,  although  Zames 
[44]  has  put  forth  some  compelling  arguments  in  favor  of  the  L* 
norm  for  sensitivity  optimization. 

However,  the  existing  L°°  norm  robustness  optimization 
results  in  [4]  and  [  40  ]  are  inadequate  for  several  reasons. 
First,  the  problem  formulations  use  only  sensitivity  or  only 
stability  margin  as  a  performance  measure,  making  no  provision 
for  the  trade-offs  between  the  two  which  are  the  salient  fea¬ 
ture  of  most  practical  feedback  designs.  Also,  the  L“-norm 
sensitivity  optimization  results  [4 0 ]  are  limited  to  single — 
loop  feedback  systems.  The  L°°-  norm  stability  margin  results 
[4]  are  only  slightly  better,  being  limited  to  decoupled  (i.e., 
diagonal)  multiloop  closed-loop  transfer  matrices.  Addition¬ 
ally,  for  both  L*-norm  sensitivity  optimization  and  for  L^-norm 
stability  margin  optimization,  the  solutions  in  [  4]  and  [  40] 
become  degenerate  and  fail  to  exist  in  the  important  case  where 
the  open-loop  plant  has  poles  or  has  zeros  on  the  imaginary 
axis.  For  example,  "type-1"  systems  (i.e.,  with  integral  feed¬ 
back)  are  beyond  the  scope  of  the  existing  stability  margin  op- 
mization  results.  Of  these  limitations,  perhaps  the  failure 
to  allow  trade-off  between  minimal  L^'-norm  sensitivity  and 
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maximal  L  -norm  stability  margin  is  the  most  serious:  virtually 
all  practical  design  problems  require  trade-offs  between  sensi¬ 
tivity  and  stability  margin.  It  is  our  expectation  that  theore¬ 
tical  links  between  this  problem  and  the  superficially  unrela¬ 
ted  but  mathematically  similar  problems  of  complex  interpolation 
and  optimal  Hankel  model  reduction  [.9, 41 , 47]  will  prove  useful 
in  this  effort. 

This  avenue  of  research  proved  to  be  spectacularly 
successful  this  year  when  we  solved  the  multivariable  L°°-norm 
sensitivity  optimization  problem  using  this  approach.  The 
results  are  briefly  summarized  in  [47],  and  a  more  detailed 
paper  is  currently  being  prepared  for  presentation  at  the 
June  1983  American  Control  Conference.  We  are  quite  pleased 
with  this  success  because  the  multivariable  problem  had 
resisted  the  efforts  of  top-notch  researchers  for  several 
years  [40,  44]. 

Perhaps  more  importantly,  our  solution  to  the  multivari¬ 
able  L®  sensitivity  optimization  problem  may  well  open  the 
door  for  the  optimal  solution  of  realistic  L®  singular  value 
shaping  problems  involving  optimal  trade-offs  between 
stability  margin  singular  values  and  sensitivity  singular 
values.  Current  research  is  focusing  on  this  important  pro¬ 
blem.  In  our  opinion  success  in  this  endeavor  would  lead  to 
a  multivariable  design  tool  that  would  have  significant 
advantages  over  the  rather  imprecise  currently  available 
techniques,  e.g.,  it  would  be  preferred  over  linear  quadratic 


optimal  control  methods  for  most  applications. 


With  as  primary  motivation  the  robustness  issues,  we  are 
looking  more  carefully  at  all  of  the  norms  relevant  to  linear 
systems,  their  interrelationships,  their  connections  with  the 
Hankel  and  the  Toeplitz  operators  associated  with  the  return  dif¬ 
ference  and  the  inverse-return  difference,  and  the  topology  the 
norms  induce  in  th  space  of  transfer  matrices.  Inspiring  our¬ 
selves  from  [46  j,  it  seems  that  there  are  several  norms  rele¬ 
vant  to  linear  systems;  except  for  the  La  norm  associated  with 
LQG  theory,  each  of  these  norms  is  the  classical  spectral  norm 
of  a  Hilbert  space  operator  associated  with  either  the  system 
impulse  response  or  the  system  transfer  matrix.  Examples  of 
such  operators  are  the  Hankel  operator  and  the  Toeplitz  operator, 
the  former  being  related  to  the  Hankel  norm  and  the  latter  being 
related  to  the  frequency  response  norm.  Actually,  as  shown  in 
[46  ],  there  are  myriads  of  relevant  operators  somewhere  between 
the  Hankel  and  the  Toeplitz  structures,  and,  consequently,  there 
are  myriads  of  norms  between  the  Hankel  norm  and  the  freouency — 
response  norm;  one  is  tempted  to  speculate  that  the  topology 
becomes  stronger  and  stronger  as  we  proceed  from  the  Hankel  norm 
In  efforts  to  provide  a  more  solid  link  between  input- 
output  and  state-space  measures  of  robustness,  we  have  recently 
re-examined  linear  quadratic  optimal  control  theory  in  the  con¬ 
text  of  Hankel-Norm  balanced  state-space  coordinates.  The 
optimal  LQG  compensator  is  well  known  to  consist  of  the  cascade 
of  the  optimal  filter 
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*  A.  -1  T  A 

x  «*  Ax  +  Bu  +  r  PC  (z-Cx) 
in  m 

and  the  optimal  control  gain 
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u  =  -  r  BP  x  . 
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p  and  P  are  the  stabilizing,  positive  definite  solutions 
m  o 

to  the  filtering  and  the  control  algebraic  Riccati  equations 

0, 
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Now,  a  "balancing"  operation 

(A,B,C)  5  (TAT'1,  TB,  CT-1)  «:  (A,  B,  C) 

is  utilized  to  simultaneously  diagonalize  the  solutions  to 
the  filtering  and  control  Riccati  equations: 
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Little  though  shows  that  the  y's  are  intrinsic  quantities , 
independent  of  the  original  state  space  realization  (A,h,C) 
The  rationale  for  "balancing"  the  plant  (A,B,C)  so  as 
to  simultaneously  diagonalize  the  solutions  to  the  Riccati 
equations  is  that  the  so-defined  u's  provide  "measures"  of 
the  "importance"  of  each  state  component  x^  of  the  plant 
G(s)  in  its  balanced  realization  (A,B,C) .  This  is  easily 
seen  using  the  error  - (x-x)- covariance  interpretation  of 
p  *  diagiu, , . . . ,u  }  and  the  cost  interpretation  of  P  * 

* 


The  theory  of  this  Riccati  balancing  is  fairly  well 
developed.  Probably  the  most  striking  result  is  the  fact 

A 

that  deleting  the  compensator  state  components  x^'s 
corresponding  to  the  lowest  y^'s  leads  to  a  stabilizing  > 
reduced-order  (LQG)  compensator. 

Within  the  scope  of  the  present  research  project,  which 
is  primarily  concerned  with  robustness  and  sensitivity 
issues  in  feedback  design,  we  have  investigated  the  robust¬ 
ness  aspects  of  the  above  described  reduced-order  compen¬ 
sator  design  scheme. 

Clearly,  y^  defines  how  much  the  (balanced)  state  com¬ 
ponent  x^  is  involved  in  the  closed-loop  LQG  behavior  of 
the  system.  (In  particular,  if  y^  is  small,  then  can  be 
removed  in  the  compensator,  yet  the  design  remains  stable.) 
But  how  is  y^  affected  by  the  robustness  requirement?  The 
robustness  properties  of  the  LQG  loop  depend  strongly  on 
the  parameters  q  and  r.  In  particular,  it  is  fairly  well 
known  that  in  some  singular  situations  the  stability  margin 
is  boosted.  We  have  utilized  the  tetrahedral  truss  of 
Draper  Lab  as  a  test  bed  to  determine  how  the  y‘s  are 

affected  by  the  robustness  requirements  on  the  LQG  loop. 

This  research  was  motivated  by  an  attempt  to  derive  some 
"rules  of  the  thumb"  to  determine  which  state  coordinates 
of  a  large  space  structure  are  the  troublemakers  and  should 
be  taken  into  account  for  compensation.  The  results  art- 
reported  in  paper  (48).  A  fundamental  result  is  that  in 
most  cases  the  balanced  coordinates  are  lined  up  with  the 


modal  coordinates.  Thus  the  y ' s  can  he  interpreted  ae  the 
importance  of  the  corresponding  vibration  modes.  We  have 
observed  that  the  y ' s  of  the  vibration  modes  do  depend  on 
the  robustness  requirements  of  the  loop.  In  general,  as 
the  stability  margin  requirement  increases,  the  y's  increase 
and  the  modes  are  more  and  more  difficult  to  discriminate 
importance  wise.  In  case  of  modes  with  overlapping  eigen- 
frequencies  the  y's  increase,  as  a  warning  that  a  resonance 
can  occur.  Finally,  in  case  of  unstably  interacting  oscillators 
(in  contrast  to  stably  interacting  oscillators 3  unstably 
interacting  oscillators  can  go  180  phase  shift  under  arbi¬ 
trarily  small  perturbations  of  the  eigenfrequencies) ,  the 
y's  are  larger  than  for  stably  interacting  oscillators,  as 
a  warning  that  instability  can  occur. 
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